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ABSTRACT 

The chlorine residuals discharged in treated sewage effluents 
to streams have caused concern because of their toxic effects 
to aquatic biota and potential health hazards through 
combinations with other substances to form carcinogens. The 
concentrations of instream chlorine residuals generally 
decrease with increasing distance from the outfall due to the 
dispersion of the effluent as well as by assimilation of the 
material by chemical, physical and biochemical instream 
processes. The report deals with the development of a 
methodology applicable to the prediction of total residual 
chlorine concentration in the mixing zones of shallow streams 
and rivers. Data collected in shallow streams in Southern 
Ontario during summer 1976 and winter 1977 have been used in 
the development of prediction techniques. The field surveys 
were primarily designed to study the magnitude and extent of 
decay of the residual chlorine in streams receiving sewage 
effluents. The decay rate coefficients of the chlorine 
residuals and other parameters of models derived from this data 
are subject to approximaticais. The study is mainly concerned 
with the development of predictive techniques and thus the 
contents of the report should be viewed in this perspective. 

In one of the streams, the chlorine residuals were found to 
disappear in less than 9^5 m below the outfall; however, there 
were significant concentration variations in the transverse 
direction. The gross first order decay rates of the chlorine 
residuals in the two successive sections of the stream below 
the outfall were found to be 0.27 per minute and 0.082 per 
minute (to base 10) . The attenuation of concentration due to 
the dispersion effects was expressed by an apparent decay rate, 
assuming that the attenuation of the mean concentration of a 
conservative material in the plume follows a first order decay 
law; the apparent decay rates were found to be about 30!6 of the 
gross decay rates under survey conditions. A two dimensional 
mathematical model was developed for the prediction of the 
distribution of instream total residual chlorine concentration 
in the mixing zone of the stream. The model includes transport 
of the material by convection and decay in the longitudinal 
direction, as well as dispersion in the lateral direction. 
Some difficulties were encountered in the calibration of the 
model due to the limitations imposed by the available data, as 
well as the need to use an average value of the kinetic rate 
coefficient for the entire study stretch as required by the 
model formulation. The predictions of the model were 
intrepreted by taking the appropriate instream processes into 
consideration, thus overcoming the limitations of calibration. 
The model was applied to predict the spatial distribution of 
total residual chlorine concentration under different flow 
conditions. The predictions of the model show that the 
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residual chlorine concentration levels increase at higher flows 
at the lateral boundary of the mixing zone (one third width) at 
some downstream transects due to increased transverse spread, 
as well as reduced decay of the material. This is contrary to 
the general belief that the instream concentrations would be 
lower at higher flows. 

Investigations in another stream indicated the persistence of 
the chlorine residuals for more than 4.8 km below the outfall 
during both summer and winter seasons. Cross-sectional 
variations in concentration were present in a stretch of 610 m 
below the outfall; however, the percentage variations were 
quite small. The first order decay rates were found to be 
0.0052 per minute and 0.012 per minute (to base 10) for the two 
successive sections below the outfall during summer; whereas 
during winter the corresponding values were 0.0095 per minute 
and 0.0029 per minute (to base 10). A plug flow model was 
found to be adequate to describe the observed concentration 
distributions in this stream; the model calibration was 
satisfactory. 

The temperature dependence of the decay rate of the residual 
chlorine was expressed by the well known van't Hoff-Arrhenius 
relationship. The theoretical value of the temperature 
correction factor was found to lie in the range 1.02 to 1.0^4 at 
pH and temperature values normally encountered in practices. 
Based on the field results, a higher range of 1.12 to l.lU was 
obtained for the temperature correction factor; however, these 
values should be treated as tentative because of limitations of 
the data. 

The models described in this report are applicable to shallow 
streams in which steady state conditions prevail. For 
application to natural streams in which the hydraulic and 
biochemical processes can vary significantly from one reach to 
another, it is necessary to formulate the model such that 
reach dependent rates can be specified; further investigations 
required to develop such a model are recommended. Considering 
the limitations, the model as presented can be applied as a 
useful tool to practical situations where the reach dependency 
of the model coefficients is not significant. 
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1. INTRODUCTION 



1 . 1 GENERAL 



Chlorine is used as a disinfecting agent in wastewater treatment 
practice in order to destroy the microorganisms of public health 
concern before the effluent is discharged to a receiving body. In 
this regard, chlorination practice has been significantly effective 
in the control of the outbreak of waterborne diseases. In order to 
achieve a satisfactory disinfection of the wastewater effluents, it 
has become a comnon practice to specify a minimum residual 
concentration level in the effluent; the present Ministry of the 
Environment practice requires a minimum concentration of 0,5 mg/L 
after a contact time of 30 minutes. In recent years, however, the 
discharge of chlorine residuals to receiving water bodies has caused 
concern because of their toxic effects to aquatic biota and 
potential health hazards. A review of the toxicity of chlorine 
residuals to aquatic life and the potential carcinogenic effects 
attributable to the chlorine residuals, as well as areas of concern 
associated with conventional disinfection practices in Ontario, has 
been presented by Ellis C1976). This report deals with the 
development of assessment techniques to evaluate the impact of 
various management alternatives in the mixing zones of shallow 
streams . 

The distribution of instream chlorine residual concentrations in 
shallow streams is governed by the combined effects of chemical 
reactions, mixing of effluent with the stream water and other 
physical and biochemical instream processes. In a portion of the 
cross-section of the stream below the outfall, the concentration of 
pollutants such as residual chlorine may be higher than the 
allowable limits at some locations in the stream; this segment of 
stream is usually termed a "mixing zone" although it may be more 
appropriate to call it a "zone of non-compliance". The segment of 
stream beyond the boundaries of this zone, wherein the water quality 
is favourable for aquatic life, is usually termed a "zone of 
passage". 

The Ministry of the Environment Water Quality Criteria (1974) for 
the protection of fish and other aquatic life states that the 
discharge to receiving waters of toxic substances in concentrations 
or combinations that may be harmful to human, animal, and/or aquatic 
life, is not permitted; however, an exception is made when mixing 
zones in the vicinity of the outfall are permitted. According to 
the criteria for zones of passage and mixing, at least two thirds of 
the total cross-sectional area of a stream should be characterized 
by a water quality which is entirely favourable to the aquatic 
community at all times. 

1.2 RESIDUAL CHLORINE FORMATION 

The dissolution of chlorine in water results in the formation of 
hypochlorous acid CHOCl) and hypochlorite ion (0C1-). The 
hypochlorous acid reacts with annonia present in the wastewater 
resulting in the formation of monochloramine (NH2CI), dichloramine 
CNHCI2) and trichloramine (NCI3), Other chlorinated compounds 



may also be formed depending upon the chemical constituents of the 
stream and wastewater. A more detailed review of the fundamentals 
of chlorination is presented by Ellis (1976) and Fair and Geyer 
(1961). 

The portion of chlorine residual concentration present in the forms 
of HOCl and OCl" is referred to as the "free residual chlorine", 
whereas the portion present as chloramines is termed "combined 
residual chlorine". The sum of free and combined forms is known as 
"total residual chlorine (TRC)". In chlorinated sewage effluents, 
the residual chlorine is present as chloramines, while free residual 
chlorine is seldom observed. The allowable concentration of 
residual chlorine for protection of aquatic life is specified in 
terms of total residual chlorine. Hence, the total residual 
chlorine parameter is used in this report. 

1.3 OBJECTIVES 

This investigation was undertaken to develop a methodology 
applicable to the assessment of various management alternatives on 
the TRC levels in the mixing zones of streams below sewage treatment 
plant (STP) outfalls and thus, aid in the application of the mixing 
zone criteria. The specific objectives of the study are as follows: 

1. To determine the relative effects of dispersion and chemical 
reaction on the reduction of instream total residual chlorine 
(TRC) concentrations. 

2. To investigate the temperature dependence of the reaction rate 
of chlorine residuals. 

3. To develop mathematical models applicable to the prediction of 
spatial distribution of TRC in streams receiving sewage effluent 



4. To investigate the effects of strearaflow variations on the TRC 
levels in the mixing zones. 



II. I N5TREAM RESIDUAL CHLORINE MODELLING 
2.1 IN STREAM PROCESSES 

2.1.1 Mixing and Dispersion 

Mixing phenomenon in natural streams tends to reduce 
non-uniformities or gradients in composition or properties of a 
material in bulk and is accomplished by the movement of materials 
between various parts of the streamwater. The movement of the 
material itself occurs by a combination of bulk flow in laminar and 
turbulent regimes, and molecular and turbulent diffusion processes; 
the resultant process affecting the transport of the material is 
consDonly known as dispersion. A review of the basic physical 
exchange processes in natural streams is presented by the author 
(1974). 

The degree of mixing and dispersion in a receiving stream is 
dependent on the hydraulic characteristics of the outfall and of the 
stream channel. In the case of a shallow, uniform stream, effluent 
released from a bank outfall may hug the shoreline for a 
considerable distance. In the case of a diffused (finite width) 
outfall, the initial dilution and further mixing and dispersion of 
effluent may be dependent on the hydraulic characteristics of the 
diffuser. The channel hydraulic characteristics such as the 
cross-sectional depth and velocity distributions, presence of pools, 
rapids and meanders, etc., will influence the mixing and 
dispersion. Excessive growths of aquatic weeds and algae tend to 
reduce mixing. Ultimately, the cross-sectional concentration 
distribution attains uniformity at some distance below the outfall. 
The zone of a stream in which there exist significant variations of 
concentration in the lateral and vertical directions, is termed a 
"mixing zone". This definition is broader than the traditional 
definition of a mixing zone based on water quality criteria (i.e. 
zone of non-compliance) and will be used in this report. 

Figure 1 shows a typical plan view of a mixing zone along with 
variations in cross-sectional width, depth, velocity and 
concentration distribution of a conservative material at various 
transects below the outfall. The average concentration of the 
material in the plume at each transect, Cav» is also shown in 
Figure 1. The concentration distribution profiles indicate that 
there is an attenuation of the peak and average concentration values 
at various transects below the outfall due to transverse dispersion 
of the material. A method of evaluating the attenuation of the 
average concentration of a conservative material by mixing in the 
plume is presented in Section 2.8. 

2.1.2 Kinetics of Chlorine Deca y 

Several investigators have studied the factors affecting the decay 
of chlorine residual concentration both in lab experiments and in 
field conditions and proposed different kinetic models to describe 
the decay process. A study by Katz (1977) on the chlorine 
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dissipation and toxicity in the presence of nitrogenous compounds, 
suggests the applicability of a first order kinetic model to 
describe the decay of chlorine residuals. However, 
Hostgaard-Jensen, Klitgaard and Pedersen (1977) felt that the decay 
of chlorine in cooling waters can be best represented by assuming 
two simultaneous reactions, one being an almost instantaneous 
reaction of chlorine with a component X in water, the other being a 
second order reaction. Laboratory investigations on chlorine 
residual decay in natural stream and reconstituted water samples 
were carried by the Water Quality Laboratory, Ministry of the 
Environment CWisz, Ellis and Innis, 1977); attempts to describe the 
observed chlorine residual decay by the zero, first and second order 
kinetic models were not successful. In this investigation, a first 
order kinetic model is used to describe the decay of total residual 
chlorine concentration by chemical reaction and other physical and 
biochemical processes (except attenuation by mixing). 

2-2 MATHEMATICAL MODEL FOR MIXING ZONES 

A simplified equation applicable to the prediction of the 
concentration of chlorine residuals in mixing zones of shallow 
streams is given by the differential equation: 

9c _ 9^c , 

in which x and y are the longitudinal and transverse co-ordinates, 
respectively; c is the concentration of chlorine residual at (x,y); 
u is the mean velocity of flow in the x direction; Cy is the 
transverse dispersion coefficient; and kj is the first order decay 
coefficient. The term on the left hand side of Eq. 1 represents 
advective transport while the two terms on the right side represent, 
in order, the lateral transport by dispersion and the decay of the 
material by chemical reaction. Equation 1 is based on the following 
assumptions: 

1. The instream processes have attained a steady state and remain 
steady with time. 

2. Transport of chlorine residual by dispersion in the 
longitudinal direction is negligible compared to the advective 
transport in that direction. 

3- The stream is shallow and wide so that the vertical 

distribution of concentration is uniform at all transects. 

4. The effects of secondary currents in the lateral and vertical 
directions are negligible. 

5. The decay of residual chlorine concentration follows first 
order kinetic law, and there are no other sources and sinks of 
residual chlorine in the stream. 



Analytical solutions of Eq. 1, subject to appropriate boundary 
conditions, are readily available in the literature for bank and 
diffused outfall sources (Carslaw and Jaeger, 1956, and Lau, 1972). 
However, for streams in which there are significant depth and 
velocity variations and meanders, the effects of secondary currents 
appear to become important. Modelling techniques applicable to such 
complex situations will be considered in future studies. 

2.3 VARIATION OF HYDRAULIC PARAMETERS WITH FLOW 

The variations in channel width, depth and velocity with flow are 
expressed by the Leopold-Maddock equations. If the values of width 
b-j, depth hi, and velocity u-\, at a flow Qi , are known, then 
b2, h2 and U2 at another flow Q2 may be calculated from the 
following equations; 

b2 = biCQ2/Qi)P ... C2) 

h2 = hi(Q2/Qi)<l ... (3) 

"2 = ui(Q2/Qi)'' ... CU) 

The values of the exponents p, q and r will have to be estimated 
using field data. 

2.4 TRANSVERSE DISPERSION COEFFICIENT 

The degree of mixing in the lateral direction by turbulence and 
cross-sectional velocity variation is expressed by the transverse 
dispersion coefficient, ey. The value of ey can be determined 
from transverse concentration distribution data. In the absence of 
such data, the value may be estimated from the relation eyraehu, 
in which a^ is a dimensionless coefficient whose value varies from 
about 0.02 to 0.0^1; and h and u represent average values of depth 
and velocity, respectively. 

Lower values of a^ are indicative of lower exchange capacities of 
streams and are applicable to slow moving sluggish channels; higher 
values of cxg indicate higher exchange capacities that are typical 
of shallow fast moving turbulent streams. 

2.5 EFFLUENT PLUME WIDTH 

In the case of bank outfalls, and diffused outfalls with a diffuser 
length much smaller than the channel width, the width of a plume, 
yp, at a transect distant, x, from the source is assumed to follow 
an empirical relation expressed by the following non-dimensional 
equation: 






X 



m 



I ' i^'d^ ^x 
Ix^j i L P- X 



(5) 



in which b^ = width of channel at a transect distant x from the 
source; xj = distance from the source to a transect where the 
plume width is equal to one half of the channel width; a = 0.5 (by 



definition). The values of m and xj will have to be evaluated 
from field data analysis. Equation 5 is similar in form to the 
empirical expressions used to describe the spread of wakes and jets 
(Hinze, 1959). 

The plume width predicted by Eq. 5 represents the average plume 
boundary of a conservative material (e.g. chlorides). The plume 
width of residual chlorine will be within the value of width 
predicted by Eq, 5. 

2.6 MIXING ZONE LENGTH 

According to the definition, the mixing zone length, x^q, is the 
downstream distance between the outfall and the transect at which the 
cross-sectional concentration distribution has attained uniformity. 
An expression for mixing length'' has been developed by Fischer 
C1967), based on a theoretical analysis of turbulent diffusion 
mechanisms. Using an analytical solution of the two dimensional 
convective diffusion equation, Ruthven (1971) obtained an expression 
for the mixing zone length in shallow streams, on the assumption 
that the variation in cross-sectional concentration distribution 
does not exceed an arbitrarily chosen limit. Both of these 
approaches have resulted in an expression of the form: 



^m = 



£/u 



04n inir ' •'• *6' 



in which £ is a characteristic length dependent on the location of 
source and position of maximum surface velocity within the 
cross-section, u» is the shear velocity and oijjj is a 
non-dimensional constant. 

For the estimation of mixing length, Fischer suggested a value of 
1.8 for the constant a qj based on laboratory test data, whereas 
Ruthven obtained a value of 1.3 based on a criterion of ^0% 
variation in the cross-sectional concentration values. The 
characteristic length X. , is usually expressed as a fraction of the 
channel width b; the ratio Z/h varies from 0.5 for a source located 
at the centre of the channel to 1.0 for bank outf alls . The shear 
velocity can be calculated from the relation uft=^hS, where 
g=acceleration due to gravity and S=slope of channel bed. 



^Mixing length is defined as the distance below a point or line 
source that corresponds to the duration of an initial convective 
period in which the movement of a tracer cloud does not follow the 
one dimensional dispersion theory. 



2.7 SPATIAL DISTRIBUTION OF TOTAL RESIDUAL CHLORINE 

The concentration distribution of TRC at any transect below an 
outfall may be obtained by solving Eq, 1 subject to appropriate 
boundary conditions. For the case of a vertical line source of a 
conservative material discharged from bank and diffused outfalls 
integrated solutions of Eq. 1 are presented by Carslaw and Jaegar 
C1956), Holley, et al (1972) and Lau (1972). Since the decay of 
chlorine residual concentration is assumed to follow the first order 
kinetic law, the solutions are to be multiplied by the term 
exp(-kdx/u), where k^ is the decay rate coefficient due to 
chemical reaction and other mechanisms, except dispersion effects. 
The analytical solutions of Eq. 1 used in this study are presented 
in the Appendix. 

In case of bank outfalls, there is an immediate lateral spread of 
the plume below the source due to such factors as protrusion of 
outfall into the river, discharge velocity of outfall, etc. The 
effluent mixing in such a case can be considered to be equivalent to 
that of a diffused outfall with a finite source width, Wg, which 
has to be estimated using the field data. 

2.8 RELATIVE EFFECTS OF MIXING AND DECAY 

If the instream residual chlorine is treated as a conservative 
material, it is possible to predict the cross-sectional distribution 
of TRC due to mixing effects only by setting kd=0. A comparison 
of the predicted distribution with and without decay would then 
permit an estimation of the relative effects of mixing and decay on 
the attenuation of TRC concentrations. 

Another method of evaluating the significance of mixing relative to 
decay by other mechanisms (such as chemical reactions) on the 
reduction of TRC is to compare the rates of decay attributable to 
each of these mechanisms. The concentration prediction technique 
described in the previous paragraph does not permit a direct 
evaluation of a first order attenuation rate coefficient due to 
mixing. Therefore, an empirical method described below has been 
developed to estimate an apparent attenuation rate coefficient due 
to mixing. 

The strearaflow in the plume Qp, that dilutes the effluent at each 
transect is estimated from the equation: 

Qp = yp Qp/bx, . , . (7) 

where Qp = total streamflow. Then the average concentration of 
the conservative material in the plume at each transect is 
determined from: 

Cav = CeQe/Qp, . . . C8) 

where Ce = concentration in effluent; and Qg = effluent flow. 



Assuming that the attenuation of TRC concentration due to dilution 
follows an exponential law, an apparent decay rate coefficient due 
to mixing, kp, can be determined from a plot of (log Cgy) vs 
distance Cor time of travel). 

The gross decay rate of TRC kt, due to mixing, chemical reaction 
and other instream processes can be determined from a semi-log plot 
of the observed average TRC concentration vs. distance Cor time of 
travel). The decay rate k^, due to chemical reactions and other 
processes Cexcept mixing) can then be determined from the relation 
kd=(l<t-kp). 

It should be noted that the attenuation of concentration 
attributable to the transverse dispersion process is expressed by 
the apparent decay rate coefficient, kp, In order to evaluate the 
first order decay rate coefficient, kj. The value of kp would 
vary with the transverse dispersion characteristics and hence is 
unsuitable for prediction of concentration distribution under 
various streamflow conditions. 

2.9 TEMPERATURE EFFECTS 

The effect of temperature on the decay rate coefficients is usually 
represented by the following form of the van't Hoff-Arrhenius 
relationship CFair and Geyer, 196I): 

^2 di '^ IRT1T2 



in which k^i and kjj are the rate coefficients at temperatures 

Ti and T2 Cexpressed in °K) respectively; R is the gas 

constant (1.99 cal per OC); and E is a constant characteristic of 

the reaction termed "activation energy". The higher the value of E, 

the slower is the reaction. In practice, the value of the term 

CE/RT^T2) tends to attain a narrow range. Because of this fact, 

Eq. 9 can be written in a modified form as follows: 

^d2 = ^d^^ ^' ^ ... CIO) 

where e=CE/RT-]T2) is a temperature correction factor. The 
values of activation energy E for reactions involving aqueous 
chlorine and chloramines at room temperature usually lie in the 
range 10,000 to 20,000 cal. In sewage disinfection practice, the 
residual chlorine is mostly in the form of chloramines. Using the 
values of E at different values of pH normally encountered in 
practice (Fair and Geyer, I961) the values of 6 at various 
temperature values were computed and are summarized in Table 1 . 
Thus, the value of 6 is seen to lie in the range 1.025 to 1.031 for 
the pH range 7.0 to 8.5 and temperature range 40C to 25°C, 
normally encountered in natural streams. For practical purposes, 
the temperature dependence of chlorine residual reaction rate 
coefficient can thus be expressed by Eq, 10 by taking 9=1.03. 



TABLE 1 : TEMPERATURE CORRECTION FACTORS 



l2 



pH E 
Units Calories ^c or oq or 



7.0 


12,000 


7.0 


12,000 


8.5 


14,000 


8.5 


14,000 


8.5 


14,000 


9.5 


20,000 


9.5 


20,000 



477 


4 


477 


1.027 


498 


25 


498 


1.025 


477 


4 


477 


1.031 


498 


25 


498 


1.029 


477 


25 


498 


1.030 


477 


4 


477 


1.045 


498 


25 


498 


1.041 



2.10 COMPUTER MODEL 

A computer model HGDISEF has been developed to predict the 
transverse distribution of concentration at various transects using 
the relationships described in Sections 2.3, 2.4, 2.5, 2.7 and 2.9. 
The program incorporates point source (bank outfall) as well as 
finite width (diffused) source models; the origin of both sources is 
placed at the outfall shore. The program can be used to predict the 
concentration distribution of a conservative or non-conservative 
material, the option being determined internally by the value of the 
decay rate coefficient. If the transverse dispersion coefficient is 
zero, the model simplifies to the plug flow case. 

The input data includes the number of transects, channel hydraulic 
and plume width parameters, length of diffused outfall, transverse 
dispersion coefficient, decay rate coefficient, temperature, 
streamflow, effluent flow rate and concentration of TRC in the 
effluent. There are options built into the program to vary any 
desired input parameter. 

The output consists of a table of lateral distances from shore and 
corresponding concentration values at each transect, and includes 
predictions from both point source as well as finite width source 
models, as well as plume width values predicted from Equation 5. 
The one third width concentration values at each transect are also 
tabulated in the output matrix. 

Further details of the computer program are provided in the Appendix 
including a description of the variable names, program listing and a 
sample input and output. 
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Ill FIELD STUDIES 

3.1 DESCRIPTION OF STUDY AREA 

Investigations were carried out in four shallow streams located in 
southern Ontario. The details of these investigations are presented 
by Wisz, Ellis and Inniss (1977). Data collected in the following 
streams are used in the analysis presented in this report: 

Boyne River below Alllston STP outfall; 
Aurora Creek below Aurora STP outfall. 

Soheraatic layouts of the streams are shown in Figure 2. The channel 
alignment of the Boyne River in the study location is straight with 
the presence of a few flat bends; the alignment of the Aurora Creek 
channel is marked by the presence of several sharp bends. The bed 
material in both streams is made up of silt, sand and boulders. 

The Boyne River flows into the Nottawasaga River; the distance 
between the Alliston STP outfall and the confluence of the two 
streams is approximately 2.9 km. The mean annual daily flow in the 
Boyne River upstream of the Alliston STP is 1.62 ra^/sec and the 
minimum daily flow is reported to be 0.11 m^/sec based on 1968-7O 
data. Effluent from a food processing industry CSalada Foods) is 
discharged into the creek at about 9^5 m below the Alliston STP 
outfall. 

Aurora Creek is a tributary of the East Holland River; the channel 
length between the Aurora STP outfall and the Holland River is about 
4.8 km. Seasonal streamflow measurements in the creek indicate that 
the flow varies from 0.10 m^/sec to about 0.50 m^/sec. A small 
unnamed creek enters Aurora Creek about 1.6 km below the Aurora STP 
outfall. 

3.2 DATA COLLECTION ASPECTS 

During the field investigations, data on channel geometry, flow 
rates, tijne of travel and effluent dispersion characteristics were 
collected. Instream residual chlorine levels at various distances 
below the outfall were measured during a preliminary run. Based on 
these data, several transects were established in the stream for 
sampling purposes. At some of the transects below the outfall. 
Gross-sectional sampling points were located based on the effluent 
dispersion characteristics in order to collect data on lateral 
variations of water quality parameters including residual chlorine 
concentrations. The details of field data at each location, 
including the number of sampling runs and parameters measured are 
presented by Wisz, Ellis and Inniss (1977). A brief summary of the 
field data is presented below. 

An intensive field survey was carried out in the Boyne River during 
August 11-12, 1975. During the survey, the streamflow in the Boyne 
River upstream of the Alliston STP was measured at O.78 m^/aec and 
the average instream temperature was 21.5°C; the average STP 
effluent flow rate was 39-6 L/sec. The average concentrations of 
total residual chlorine and chloride at various sampling points for 
the duration of the survey are summarized in Table 2. 
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FIGURE 2. SCHEMATIC LAYOUT OF STUDY AREAS. 



TABLE 2: SUMMARY OF BOYNE RIVER FIELD SURVEY DATA - 
AUGUST 11-12, 1976 



Transect 
and 


Distance below 


Transverse Distance 


Avg^ 


Avg 

Chlorides ^ 


Sampling 


Outfall 


below Outfall Bank 


TRC 


as Cl- 


Station 


m 


m 


ug/L 


mg/L 


Control 






Q 


13.7 


Effluent 







1170 


131.0 


Transect A 


21 


8.Qi\^ 


100° 




Ai 




0.381 


220 


47.0 


A1 




0.76 


18 


27.0 


A2 




1.52 





16.2 


A3 




2.29 





14. U 


AM 




n.&Q 





13.6 


A5 




6.71 





13.6 


Transect B 


61 


10. 70^' 


16C 




31 




1.52 


23 


24.3 


B2 




3.05 


10 


20.4 


B3 




4.57 


2 


17.0 


BH 




7.62 





13.6 


B5 




9.^5 





13.4 


Transect C 


183 


6.40b 


3° 




CI 




1.22 


6 


20.0 


C2 




2.HH 


3 


19.0 


C3 




3.66 


1.4 


17.0 


C4 




U.88 


0.6 


15.0 


Salada Foods 


9^5 







21.0 


Outfall 










Transect D 


1148 







18.0 


Transect E 


2897 







17.0 



^Average concentration of five samples taken over a duration of 27 
hours (except three samples at Aj). 

^Channel width at the transect. 

°Mean concentration in the effluent plume. 



Intensive field investigations in the Aurora Creek were carried out 
during summer and winter. During the summer study (August 18-19, 
1976), the flow in the creek upstream of the Aurora STP outfall was 
0.12 m3/3ec and the average instream temperature was 18.3^C. 
The average effluent flow rate was 98 L/sec. The flow in the 
unnamed creek was not measured; however, based on a mass balance of 
conductivity, the flow in this creek was estimated as 0.017 m3/3ec 
whereas a mass balance of chlorides resulted in an estimation of 
0.061 m3/3ec. During the winter study (February 9, 1977), the 
flow in the Aurora Creek was 0.166 m3/sec including an effluent 
flow of 87.5 L/sec; the average instream temperature was 5.5*^C. 
There was some ice cover in the lower sections of the creek. The 
average concentrations of TRC, conductivity and chlorides at various 
transects are summarized in Table 3 for the summer studies. The 
cross-sectional mean concentrations of TRC obtained during the 
winter study are also presented in Table 3- 
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TABLE 3: SUMMARY OF AURORA CREEK FIELD SURVEY DATA 



Transect and 


Distance below 




Summerf August 


18- 


-19, 1976 Winter: February 9, 


1977 


Sampling 


Outfall 


Channel Width 


TRC 


Conduct iv 


Lty 


Chloride 


TRC 




Station 


m 


m 


ug/L 


umho/cm 




mg/L 


ug/L 




Control 









626 




ff 







Effluent 







1536 


2190 




342 


2330 




ft 


15.2 


1.83 


722 


1353 




187 


1379 




ii 


91. U 


3-50° 


615 


1362 




169 


918 




#- 


U02.3 




520 


1434 




172 


980 




© 


643.7 




404 


1380 




208 


902 




i. 


1127.0 




287 


1400 




180 


824 




r 


1609.0 




180 


1415 




1ft 


786 




§ 


2092.0 




82 


1354 




tfis 


607 




t 


2897.0 




11 


1332 




173 


403 




J 


4828.0 




0.8 


1270 




160 


213 





^Average concentration of seven samples taken over a duration of 19 hours. 
^Average concentration of four samples taken over a duration of 10 hours. 
•^Estimated average channel width at each transect. 



IV DATA ANALYSIS AND MODELLING RESULTS 

U.l BOYNE RIVER BELOW ALLISTON SIP OUTFALL 

To investigate the effect of streamflow variations on the total 
residual chlorine (TRC) distribution, the following streamflow 
values ranging from 0.25 m3/sec to 1.64 m3/sec were used. The 
width, depth and velocity values at various flows were computed 
using the following values of the exponents in Equations 2 to H: 
p=0.1; q=0.5 and r=0.J4. The flow and associated width, depth and 
velocity values are presented in Table 4. Computation of the 
transverse dispersion coefficient is based on ^^=0.02. The 
predicted mixing zone length values at the above mentioned flows 
with c(jjj=i.8 and Manning's roughness coefficient n=0.035, are 
presented in Table 4 for ^= b/2, J^=2b/3 and ^=3b/4. The mixing zone 
length is found to decrease with an increase in the streamflow, 
reflecting the effect of increase in transverse diffusion at higher 
flows. 

In order to determine the gross TRC decay rate coefficient k^ (due 
to mixing and all other mechanisms), a semi-log plot of observed 
average TRC concentration vs. distance was prepared (see Fig. 3). 
An examination of the plot indicates that the rate coefficient 
varies significantly in the two sections as indicated by the slopes 
of the solid lines. The estimated rate coefficients (kt) for the 
two sections are 0.27 and 0.082 per minute (to base 10), 
respectively, at the observed mean temperature of 21.5*^0. 

The plume width values were determined using the transverse 
distribution of chloride concentration at various transects; the 
values of xj and ra were evaluated to be 59.4 m and 0.43 
respectively, from a log-log plot of (yp/b^) vs. x. 

By treating the average TRC as a conservative material, the average 
TPC concentration values in the plume due to mixing effects only, 
computed from Eqs. ? and 8, are: 177.69, 112.65, 70.38 and 56,98 
ug/L at the transects A, B, C and D respectively. Using these 
values, the apparent attenuation rate coefficients due to mixing 
(kp) are found to be 0.069 and 0.023 per minute (to base 10) for 
sections AB and BC, respectively (see Fi^re 3), representing 
approximately 25% and 28$ of the gross decay rate values. Thus, 
attenuation of TRC concentration due to mixing only accounts for 
about 30$ of the observed decay. The corresponding decay rates of 
TRC (k^) attributable to chemical reactions and other mechanisms, 
determined from k(j=(kt-kp), are 0.20 and 0.059 per minute (to 
base 10). Since a different decay rate for each section could not 
be specified for the model used, an average value of kjj (obtained 
from model calibration) has been used in the computations. 

In order to calibrate the model, several modelling runs were carried 
out with various combinations of kj »^e ^"^ "s* Based on 
comparisons between the observed values of TRC and predictions of 
the model at various transects, it was found satisfactory to use 

0(e=0.02 and W3=1.2 m. The modelling results also showed that the 
predicted values of both models were higher than the observations at 
Transect A; this may be attributed to the effect of initial mixing 
mechanisms since Transect A is too close to the source, being 
located at 21.3 m below the outfall. Therefore, the observed TRC 
data at Transect A were felt to be unsuitable for calibration. The 
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TABLE H: BOYNE RIVER - HYDRAULIC PARAMETERS AND MIXING ZONE LENGTHS AT VARIOUS FLOWS 



Flow 



Width 
m 



Depth 

ID 



Velocity Mixing Zone Length", m 
m/sec i=h/2 £=2b/3 £ =3b/^ 



Remarks 



0.20 


8.60 


0.183 


0.131 


1249 


2220 


2813 




0.41 


9.27 


0.259 


0.172 


1086 


1931 


2445 




0.82 


9.97 


0.366 


0.227 


925 


1646 


2084 


. . . Summer Survey 
Conditions 


1.64 


10.59 


0.518 


0.300 


796 


1415 


1791 





^Using Manning's Roughness Coefficient, n=0.035. 
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FIGURE 3. ATT|fNUATION AND DECAY OF TOTAL RESIDUAL CHLORINE 
IN BOYNE RIVER 
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selection of a value of kj posed some problems. The decay rates 
vary significantly in two successive reaches (see Fig. 3). However, 
the formulation of the model assumes constant values of k^ and 
Sy in all the reaches, necessitating the use of a value for each 
coefficient that is representative of the entire study stretch of 
the stream. When a higher value of k^ (0.20 per minute) was used, 
there was a better agreement between the observations and 
predictions at Transect B, while the predictions at Transect C were 
too low. After several trials with different values of k^ (within 
the range of values stated in the previous paragraph), reasonable 
agreement between the observed and predicted values was obtained at 
either Transect B or Transect C, but not at both transects 
simultaneously. Based on these analyses, kd=0-12 per minute (to 
base 10), was selected to be a representative value for the entire 
study stretch. With the combination of the parameters ae=0.02, 
W3=1.2 m, and kcj=0.12 per minute, the TRC values predicted by 
both point source and finite width source models were found to be 
generally higher than the observed values at Transects A and B, 
whereas the predicted values at Transect C were in reasonable 
agreement with the observations. For the purpose of further 
analyses, the finite width source model predictions were used; the 
observed and predicted values of TRC are shown in Figure UCc). 

The difficulties encountered in the calibration procedure are 
attributable to the assumption of constant values of kj and ey 
in the model formulation. In shallow streams with highly variable 
morphological characteristics, the hydraulic and biochemical 
processes are subject to rapid changes from one reach to another, 
thus rendering the assumptions of constant kj and Cy 
inappropriate. This analysis brings out the need for the 
development of a model which allows for such variables. In spite of 
the limitations outlined above, the model can still be applied as a 
useful tool to many practical situations with the exception of those 
cases where the variabilities are significant, e.g. initial stream 
reaches immediately below the outfalls. 

The TRC values predicted by the finite width source model for 
various streamflow values are presented in Figures iJ(a)-(d). From 
these plots, it can be seen that the peak concentration of TRC 
decreases with an increase in the streamflow. The plots also 
indicate that the width of plume at each transect increases with the 
streamflow reflecting the effect of increase in transverse mixing 
rates at higher flows. 

The variation of TRC value at the one third width point of each 
transect (from the outfall shore of channel) with streamflow and 
downstream distance, shown in Figure 5(a), indicates that the TRC 
value increases with an increase in the flow, even though the 
dilution ratios of streamflow to SIP effluent flow are higher. 
However, as the streamflow tends to zero, the TRC values at 
transects close to the outfall can be expected to increase. The one 
third width TRC values are seen to have highest values at Transect 
C. The channel is narrow at this transect, having a width of 6.M m 
(Table 2). If the channel width is increased to 10 m at a 
streamflow of 0.82 m^/sec, the model predictions result in 
considerably lower values of the one third width TRC values as shown 
in Figure 5(b). Thus, variations in channel width, which are quite 
coninon in natural streams, appear to have an effect on the 
transverse distribution of the TRC concentrations. 
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FIGURE 5(b) EFFECT OF INCREASED CHANNEL WIDTH AT TRANSECT 
ON ONE THIRD WIDTH PREDICTED TOTAL RESIDUAL 
CHLORINE VALUES IN BOYNE RIVER. 



4.2 AURORA CREEK 3EL0W AURORA STP OUTFALL 
4.2.1 Summer Study 

Total residual chlorine concentrations were predicted under four 
streamflow conditions. The flow rates and related hydraulic 
characteristics of the channel are summarized in Table 5. (In 
Equations 2 to 4, the values of exponents used are p=0.1, q=0.5 and 
r=O.U). In the plume width computations, the values of xi and m 
were taken equal to 30.5 m and 0.9 respectively; this is based on 
the fact that the effluent (discharged from the bank outfall) 
spreads across the entire channel width in a distance of about 60 m 
below the outfall. The significance of transverse mixing was 
investigated by taking ag=0.02, as well as ctg^O.O (plug flow 
model). The results indicated that the effluent plume spreads 
across and fills the entire channel width at about 90 m below the 
outfall, and the cross-sectional variations in concentration are 
within 10$ at about 500 m below the outfall. The computed mixing 
zone length values using cc^,= i .8, Manning's roughness coefficient 
n=0.025 and i=b/2, i= 2b/3 as well as £=3b/4, are presented in Table 
5. The mixing zone length generally decreases with an increase in 
streamflow, as expected. 

A semi-log plot of the cross-sectional mean of observed values of 
TRC vs. distance below the outfall is shown in Figure 6. The rate 
of change of concentration downstream of the outfall is seen to 
change significantly below Transect F. The decay rate coefficient 
is found to be 0.0052 per minute (to base 10) for the reach AF and 
0.012 per minute (to base 10) for the reach FH. The decay rate 
below Transect F may have been affected by the inflow from the small 
creek. For modelling purposes, a decay rate value of 0.0052 per 
minute (to base 10) has been used. No attempt was made to evaluate 
the apparent attenuation rates of TRC due to mixing since the plume 
spreads across the entire channel width in a relatively short 
distance below the outfall. 

The model calibration procedure resulted in TRC values predicted by 

the diffused source model being too low (for WgzSl cm and 
Cte=0.02); whereas, point source model predictions were in good 

agreement with the observations at Transects B to F. (Note: When 
ae=0. diffused source model predictions will be the same as those 

of the plug flow model). 

As mentioned earlier, point source model predictions are based on 
o.Q=:Q.02, as well as Q. q=0 (plug flow model). The predictions 
indicated that at Transects B to E, transverse variation of TRC 
concentration could exist, although the effluent has filled the 
entire cross-section. The longitudinal distribution of TRC 
predicted by the plug flow model for the survey condition presented 
in Fig. 7 shows good agreement with the observed values at Transects 
B to F. The difference between the predicted and observed TRC 
values at Transects G and H may be due to the effects of inflow from 
the small creek on the decay rate coefficient. Plug flow model 
predictions at other flows are also presented in Fig. 7. No attempt 
was made to predict the one third width TRC values under different 
flow conditions since the cross-sectional concentration variations 
become too small in a relatively short distance below the outfall. 
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TABLE 5: AURORA CREEK - HYDRAULIC PARAMETERS AND MIXING ZONE LENGTHS AT VARIOUS FLOWS 



Flow 
m3/3ec 



Width 
m 



Depth 
m 



Velocity 
m/sec 





0.10 


3.23 


0.173 


0.177 




0.17 


3.41 


0.226 


0.216 




0.22 


3.51 


0.259 


0.2UU 




0.3M 


3.66 


0.321 


0.290 




O.iiS 


3.76 


0.371 


0.325 


^ 











Mixing Zone Length*, m 
£=tV2 p=2b/3 i?.=3b/4 



•Using Manning's Roughness Coefficient, n=0.035. 
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FIGURE 6. DECAY OF TOTAL RESIDUAL CHLORINE IN AURORA CREEK. 
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4.2.2 Winter Study 

The possible existence of cross-sectional concentration variations 
between Transects A to D Cas predicted by the model) was confirmed 
by the data collected during this study. The decay rate coefficient 
vBs evaluated from a semi-log plot of cross-sectional mean value of 
TRC vs. distance below the outfall (Fig. 6). The rate of change of 
concentration is seen to increase significantly below Transect F as 
well as below Transect H, as in the case of the summer study. 
However, the concentration values at Transects A and B deviate from 
the line of best fit. The decay rate is 0.00095 per minute for the 
reach AF and 0.0029 per minute for reach FH (to base 10). 
Prediction of the instream TRC levels under various flow conditions 
was not undertaken. 

4.2.3 Temperature Effects 

The value of the temperature correction factor was evaluated from 
Eq. 10 using the results of the sunmer and winter studies. The 
values of 6 were found to be 1.14 and 1.12 based on the decay rates 
of the reach AF and of the reach FH, respectively. These values are 
much higher than the theoretical values given in Table 1 . However, 
these values should be treated as tenative since they are derived 
from limited data. 

4.3 DISCUSSION OF RESULTS 

The evaluation of model parameters and calibration of the model are 
based on data of one field survey at each location. The value of 
transverse dispersion coefficient was assumed because of limited 
data on effluent dispersion characteristics. The calibration of the 
model posed some problems due to the necessity to use the same 
values of some of the coefficients for all reaches to conform with 
the assumptions involved in the formulation of the model, as well as 
some of the limitations associated with the available data. In case 
of the Boyne River, the model predictions satisfactorily represented 
the observed values of TRC at Transect C only; there were 
significant differences between the predicted and observed values of 
TRC at Transects A and B. In case of the Aurora Creek, there was 
good agreement between the observed and predicted values of TRC at 
Transects B to F; the model predictions did not satisfactorily 
represent the observed data at Transects A, G and H. The effects of 
the characteristics of stream and wastewater, stream turbulence, 
etc. on the decay of instream chlorine residuals also need to be 
investigated and included in the modelling procedure. These 
limitations should be taken into consideration in the interpretation 
and use of the results presented in this report. Nevertheless, this 
investigation has unveiled some interesting results which have 
significant implicatic^s on the management of water quality in 
effluent mixing zones. A discussion of some of the important 
aspects is presented below. 

The expression for effluent plume width (Eq. 5) is written in 
non-dimensional form to account for possible variations in the plume 
width in proportion to the channel width in natural streams. Field 



27 



data collected during effluent dispersion studies in the Grand River 
below Waterloo ,STP indicate that there is some justification for the 
above-mentioned assumption. However, a modification of the 
expression by including transverse dispersion parameters, and 
further testing and validation may be required. 

Using the concentration prediction results with k^=0 (conservative 
material), it is found that the cross-sectional concentration 
distribution approaches uniformity at about 915 m in the Boyne River 
and 500 m in the Aurora Creek. These estimates are found to agree 
well with the predicted mixing zone length values of 925 m with 
Jtb/2 in the Boyne River, and 490 m with £=3b/4 in Aurora Creek. 
Thus the ratio H/h in each stream has a different value although 
both streams have bank outfalls. 

The effect of transverse mixing in the Boyne River is estimated to 
cause about 30it diminution of average plume concentration of TRC 
under survey conditions. Since the observed concentration 
distributions result from the combined effects of dispersion and 
chemical reaction processes, it is essential that the effects of 
dispersion on the concentration values are separated when kinetic 
models are used to evaluate the decay rate coefficients using the 
observed data. 

The finite width model predictions of instream TRC concentrations in 
the Boyne River indicate that at higher flows, residual chlorine is 
propagated to greater downstream distances. Referring to Figure 5, 
the one third width TRC value is seen to increase at some transects 
at higher flows; this is contrary to the belief that at higher flows 
dilution ratios are higher and hence, concentrations are lower. 
This phenomenon can be attributed to the following facts: At a 
higher flow the time of travel to a given transect is lower; 
consequently, the mass of residual chlorine remaining at that 
transect is higher since the decay of mass due to chemical reaction 
is proportional to exp (-k^jt) . However, the transverse 
coefficient is higher at a higher flow, resulting in an increased 
transverse spread of the material present in the stream. Thus, the 
one third TRC value increases at a higher flow due to combined 
effects of increased mass remaining and increased transverse 
spreading. Ultimately, the one third width TRC value would start to 
decrease at some downstream point at which reduction of the mass by 
chemical decay becomes significant. This analysis suggests that 
higher strearaflows do not necessarily provide safer zones for fish 
passage. 

The variations in channel width at various transects may also affect 
the lateral distributions of concentration, and hence, the lateral 
boundary of the zone of passage. In the Boyne River, the channel is 
narrower at Transect C compared to the widths at other transects 
(Table 2). Referring to Figures 5(a) and (b), it can be seen that 
the one third width predicted TRC values would be significantly 
lower at Transect C if the channel at the transect is widened to 
10 m from 6.4 m at a streamflow of 0.82 m^/sec. There are several 
assumptions inherent in an analysis of this nature, and hence the 
results have limitations. However, it is presented in this report 
to identify the importance of an aspect of natural stream morphology 
in water quality investigations. 
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In the Aurora Creek, the effluent plume was found to spread across 
the entire channel in a short distance below the outfall during the 
summer survey; the cross-sectional concentration distribution is 
likely to have attained uniformity in a distance of 500 m below the 
outfall (based on mixing zone length values). Thus, it is 
reasonable to assume that the one third width TRC values are very 
close to the cross-sectional average values at all transects 
downstream of the predicted mixing zone length. Referring to Figure 
7, the predicted average TRC concentration values are seen to 
decrease with increased flow up to Transect H, i.e., 2897 m below 
the outfall. 
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V SUMMARY AND CONCLUSIONS 

The report outlines the development of a methodology applicable to 
investigate the persistence of total residual chlorine (TRC) 
concentrations in streams receiving sewage effluent, including two 
case studies. The relative effects of dispersion and chemical 
reaction processes on the reduction of instream TRC values are 
examined. The application of two dimensional steady state models to 
predict instream TRC levels is described. The effect of variation 
of streamflow on the TRC levels at the lateral boundary of the 
mixing zone has been investigated. The temperature dependence of 
decay rates of TRC has been examined. The assumptions and 
limitations of the techniques involved in the methodology are 
discussed. 

A listing of the computer program HGDISEF used in this investigation 
is included in the Appendix. The program includes both point source 
(bank outfall) and finite width source (diffused outfall) models. 
The data requirements and typical input and output of the model are 
also included In the Appendix. 

The field data collection methodology was generally satisfactory; 
however, it would be desirable to collect data at more transects for 
more detailed modelling purposes. Further details on equipment and 
data collection aspects are described by Wisz, Ellis and Innis 

(1977). 

Using the methodology presented in this report, it is possible to 
predict the magnitude and extent of persistence of instream residual 
chlorine concentration at various streamflow and temperature 
conditions. The implications of various alternatives available for 
reducing TRC concentrations can be evaluated using the predictive 
techniques presented here. Such procedures would aid the selection 
of a suitable method of minimizing the impact of the residual 
chlorine on the water quality of receiving streams. The models 
presented in this report are useful tools in the evaluation of 
various options involved in the design of new outfalls in streams 
and rivers. 

The following conclusions are made based on the results of this 
investigation. 

1. The chlorine residual in the Aurora Creek is found to persist 
for about U.8 km below the outfall during summer as well as 
winter; whereas in the Boyne River, the residuals disappeared 
in less than 9^5 m below the outfall during summer. 

2. Transverse dispersion of effluent in the Boyne River is 
estimated to cause a reduction of about 30? of the average TRC 
concentration in the effluent plume under survey conditions. 
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3- The observed first order decay rates in the Boyne River during 
sumner were 0.27 and 0.082 per minute (to base 10) in two 
successive reaches; in the Aurora Creek, the decay rates were 
0.0052 and 0.012 per minute (to base 10) during the summer, and 
0.0095 and 0.0029 per minute (to base 10) during the winter, in 
two successive reaches of the stream. The decay rates in the 
Boyne River were obtained from limited data and hence, they 
should be considered tentative. 

k. Model predictions indicate that at some downstream transects, 
the TRC concentration at the lateral boundary of the allowable 
zone of non-compliance (i.e. one-third width) could increase at 
higher flows, contradicting the belief that the concentrations 
would be lower at higher flows. This is attributed to the 
decrease in time of travel and increase in transverse 
dispersion at higher streamflow. 

5. The temperature correction factor 6 for TRC, determined from 
the van't Hoff-Arrhenius relationship, is found to vary from 
1.025 to 1.0^5 in a temperature range of U to 250C and a pH 
range of 7 to 9-5 units. Based on field data, the temperature 
correction factor is estimated to be 1.12 to 1.14; however, 
these values may not be representative due to limitations in 
the field data. 

6^ The model needs to be improved by including reach dependent 

dispersion and reaction rate coefficients in order to overcome 
the difficulties encountered during the model calibration. 
However, the present model can still be applied to streams in 
which the coefficients in different reaches do not vary 
significantly. 
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VI RECOMMENDATIONS 

The results of this investigation have shown that dispersion 
characteristics have a significant effect on the attenuation of 
instream TRC levels. The results of field investigations and data 
analysis indicate that the decay rate coefficients and the 
persistence of instream TRC concentrations vary considerably from 
one stream to another. This may be attributed to the 
characteristics of sewage effluent and of stream water, and instream 
processes including dispersion, reaeration and the reaction of 
chlorine with such compounds as ammonia, phenol, etc. Temperature 
is one of the most important factors affecting the reduction of 
instream TRC concentrations. The distribution of instream TRC 
concentration is also influenced by the morphological and hydraulic 
characteristics of the channel including variations in streamflow 
caused by withdrawal of water, tributary inflow, ground water 
recharge and discharge, etc. In order to develop predictive models 
which include all of these processes and factors, it is recommended 
that the following investigations be carried out: 

1. Investigation of the effects of sewage effluent and stream 
water quality characteristics on the reduction of TRC levels. 

2. Investigation of the effect of hydraulic characteristics of 
streams on the kinetic rate coefficients. 

3. Further studies on the temperature dependence of the decay rate 
coefficients. 

4. Development of predictive models applicable to natural streams 
to take into account the results of the above-mentioned 
investigations, and morphological and hydraulic characteristics 
of the channel. 
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APPENDIX 

A. Analytical Solutions 

B. Computer Program - HDGISEF 

B.I Description of Variable Names 

B.2 Units 

B.3 Listing of Program 

B.U Sample Input and Output 
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A. ANALYTICAL SOLUTIONS 
Solution for Point Source ' 



c{x,y) = exp^ 

bhu \ u 



CO 



- I-2E 



COS( JTry/b) 



7f (l + <^j) 



1/2 



exp 



ll-{,-(i.*,)-} 



where (/>■ = 

^ bu 



2J7rey ^'^ 



Solution for Finite Width Source^ 



U1 



c(x,y) = ;rr~~. expl Kerf 



r -d^ \ I . /Ws-y\ ./Ws + y 



CO 



2 h u Ws \ u 



^ 



Kerf -V-WE 



'Z' 



n=l 



erf 



'2nb + Ws+y \ /2nb-W3+y 



^ 



-erf 



^ 



/2nb+Ws-y\ /2nb-Ws-y 

+ erf ^ I - erf 



^ 



V—)^ 



where ^ - 2^xey/u 



and erf I lis the error function defined by 



erf(0 = ^ fe-''^6y 



The following notations are used in the analytical solutions 

b=channel width 

o=conoentration in stream at a point (x,y) 

C =concentration in the effluent 
e 

e =transver3e dispersion coefficient 

y 

h=depth of water in stream 

k.=kinetic rate coefficient 
a 

n= number of images 

Q =effluent flow rate 
e 

u=velocity of streamflow 

w rwidth of diffused outfall 
s 

x= longitudinal distance below outfall 
y=lateral distance from outfall bank side 
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B. COMPUTER PROGRAM - HGDISEF 

B.I DESCRIPTION OF VARIABLE NAMES 

Input Variables : 

TITLE: Title of study 

NTR: No. of transects 

NY: No of strips of equal width at transects 

X: Longitudinal distance to a transect below outfall, x 

BST: Channel width at a transect at reference streamflow 

QRS: Reference streamflow 

XHALF: Distance below outfall at which plume width is one-half of 

channel width, x-, 
PLN: Exponent in plume width equation, m 
BS: Average channel width of study stretch at reference 

streamflow, b 
HS: Average channel depth of study stretch at reference 

streamflow, h 
US: Average velocity at reference streamflow, u 



BEX 
HEX 
UEX 



Width exponent, p 
Depth exponent, q 
Velocity exponent, r 



WS: Length of diffused outfall, w 

COEFEY: Non-dimensional coefficient of lateral dispersion 

equation, a 
s 

COEFT: Temperature correction coefficient, 6 
TEMPS: Temperature, C at which reaction rate coefficient is 
known 



RKS: First order reaction rate coefficient, k. 

QRVR: Design streamflow 

TMP: Design temperature 

CSPNT: Effluent concentration, C 

' e 

QSTP: Effluent flowrate, Q 
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Output Variables : 



NO: Serial No. (Col. 1) 

Y: Transverse distance from shore (Col. 2) 



CPT 
RPT 
CDF 
RDF 

EY: 



Predicted concentration - point source model (Col. 3) 
Patio of CPT to maximum concentration at the transect (Col. 4) 
Predicted concentration - finite width source model (Col. 5) 
Ratio of CDF to maximum concentration at the transect (Col. 6) 

Transverse dispersion coefficient 
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B.2 UNITS 

Length) 

Width ) -meter 

Depth ) 

Velocity - meter per second 



3 
Streamflow - ra per second 



STP flow - m^ per second 

Transverse dispersion coefficient - m per second 

Reaction Rate - per second, to base e 



Temperature - C 



Concentration - mg/1 or ug/1 or other units 

Note: The unit of output concentration is the same as the unit of 
input concentration 



39 



B.3 LISTING OF PROGRAM 
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B.i» SAMPLE INPUT AND OUTPUT 



INPUT 

liCr-'TtH:: OR^^j /■HHLK-KU^ bir- MS? iJt? tCX :• Ht:-- k LtX 
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